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We study random lattice networks consisting of resistor like and diode like bonds. For investigat- 
ing the transport properties of these random resistor diode networks we introduce a field theoretic 
Hamiltonian amenable to renormalization group analysis. We focus on the average two-port re- 
sistance at the transition from the nonpercolating to the directed percolating phase and calculate 
the corresponding resistance exponent <j> to two-loop order. Moreover, we determine the backbone 
dimension Db of directed percolation clusters to two-loop order. We obtain a scaling relation for 
Db that is in agreement with well known scaling arguments. 
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Percolation is a leading paradigm for disorder. It pro- 
vides an intuitively appealing and transparent model of 
the irregular geometry which occurs in disordered sys- 
tems. Moreover, it is a prototype of a phase transi- 
tion. Though the usual isotropic percolation (IP) has 
attracted most attention, directed percolation (DP)J^] 
is a sexy model for study as well. DP shows many 
qualitatively new features not appearing in IP. DP is 
perhaps the simplest model resulting in branching self- 
affine objects. It has many potential applications, in- 
cluding fluid flow through porous media under grav- 
ity, hopping conductivity in a strong electric fieldjU, 
crack propagation Q, and the propagation of surfaces 
at depinning transitions]^]. Furthermore, it is related 
to epidemics with a bias[^| and self-organized critical 
models Q. While the transport properties of IP have 
been studied extensively (|, ||, [l(], [ll], [l^, relatively 
little is known about transport in DP. The transport 
properties of DP have not been addressed hitherto by 
using sophisticated analytic methods like renormalized 
field theory. 

A model which captures both, IP and DP, is the 
random resistor diode network (RDN) introduced by 
Redner |]l4|, [l5| , [To]] . A RDN consists of a (/-dimensional 
hypercubic lattice in which nearest-neighbor sites are 
connected by a resistor, a positive diode (conducting 
only in a distinguished direction), a negative diode (con- 
ducting only opposite to the distinguished direction), 
or an insulator with respective probabilities p, p+, p_, 
and q = 1 — p — p+ — p—. In the three dimensional 
phase diagram (pictured as a tetrahedron spanned by 
the four probabilities) one finds a nonpercolating and 
three percolating phases. The percolating phases are 
isotropic, positively directed, or negatively directed. Be- 
tween the phases there are surfaces of continuous tran- 
sitions. All four phases meet along a multicritical line, 
where < r := p+ = p_ < 1/2 and p — p c { r )- On 
the entire multicritical line, i.e., independently of r, one 
finds the scaling properties of usual isotropic percola- 



tion (r= 0). For the crossover from IP to DP see, e.g., 
Hef.0. 

In this letter we focus on the vicinity of the critical 
surface separating the nonpercolating and the positively 
directed phase. Here, typical clusters are anisotropic 
and they are characterized by two different correlation 
lengths: £|| (parallel to the distinguished direction) and 
£_l (perpendicular to it). As one approaches the criti- 
cal surface, the two correlation lengths diverge with the 
exponents and v\_ of the DP universality class. 

In the first part of this letter we study the average 
resistance between two connected sites x and x' when 
an external current / is injected at x and withdrawn at 
x'. We choose n = l/vd(l, . . . , 1) for the distinguished 
direction. We assume that the bonds t>uj\ between two 
nearest neighboring sites i and j are directed so that 
bu j\ ■ n > 0. The directed bonds obey the non-linear 
Ohm's law 



o-ij {V S - Vi) [Vj - Vi] = h 



(1) 



where Vi is the potential at site i and Iij denotes the 
current flowing from j to i. The bond conductances 
tJij are random variables taking on the values a, a0 (V), 
cr0(—V), and with respective probabilities p, p+, p-, 
and q. a is a positive constant and 9 denotes the Heav- 
iside function. Note that the diodes are idealized: un- 
der forward-bias voltage they behave as "ohmic" resistors 
whereas they are insulating under backward-bias voltage. 
A central role in our theory is played by the power 



P ({V}) - £ (V 3 - Vi) [V 3 ~ vx 



(2) 



<«,J> 



dissipated on the network. The sum in Eq. (Q) is taken 
over all bonds of the lattice. Following an idea by 
Stephen |l8[ and its generalization to networks of nonlin- 
ear resistors by Harris |lj| we exploit correlation functions 
of ip\ (x) = exp (i\V x ) as generating functions of the re- 
sistance R(x,x') between x and x'. Note that A = il 



2 



is an imaginary current. With help of the saddle point 
method (the integration is not Gaussian due to the 
functions) we find 



(1> X (x)V-A (*')> 
1 

~ z 

X 2 



ex exp 



I J <M ', exp 

R (x, x) 



-^P({V})+iX(V x -V x ,) 



(3) 



provided that the condition I 2 3> o holds. Z in Eq. (||) 
stands for the usual normalization. 

We are interested in the average (•••)c of R over 
all diluted lattice configurations which we will denote 
by Mr. Hence we switch to D-fold replicated voltages 



Vi 



■■,Vi 
> \ 



(D. 



and imaginary currents 
The replication procedure 



induces the effective Hamiltonian 



In yexp 



-P 



(4) 



For technical reasons [g0| we switch to discretized volt- 
ages 9 and currents A taking values on a discrete D- 
dimensional torus. For the saddle point method to be 
reliable we work near the limit when all the components 
of A are equal and continue to large imaginary values. Ac- 
cordingly we se t[jl9t A^"' = iAo + ^ a - > with real positive 
Ao and ^ a ', Yl a =i = ^> an< ^ impose the conditions 
A 2 , < D- 1 and ^ < 1. 

To refine H rcp towards a field theoretic Hamiltonian, 
we expand H rcp in terms of t/>? (a;) . The steps are analo- 
gous to those in Ref.Jl^] and are skipped here for brief- 
ness. The so obtained expression is converted into a 
Landau-Ginzburg- Wilson-type functional 



H 



ddx oE^w 



r - V 2 + wX 2 



A^O 



+ (fl(Ao)-fl(-Ao))v-V]^(x) 

A*,A*',X+A*'#0 

by applying the usual coarse graining procedure. The pa- 
rameter r specifies the "distance" from the critical sur- 
face under consideration. The vector v lies in the dis- 
tinguished direction, v = vn. t and v depend on the 
three probabilities p, p + , and p_. w is the coarse grained 
analog of o~ x . In the limit w — > our Hamiltonian 7i 
describes the usual purely geometric DP. Indeed TL leads 
for w — > to exactly the same perturbation series as 
obtained in[glj, ||, |§. 

We proceed with standard methods of field theory J24| 
and perform a diagrammatic perturbation calculation up 
to two-loop order. As in our previous work on trans- 
port in IP|| [l(], O, [l2|, the principle propagator 



consists of an conducting and an insulating part. Hence, 
the principle Feynman diagrams can be decomposed into 
conducting diagrams consisting of conducting and insu- 
lating propagators. These conducting diagrams can be 
interpreted as being directed networks themselves. This 
real-world interpretation leads to a substantial simpli- 
fication of the actual calculation. Instead of carrying 
out tedious summations over loop currents, we just have 
to determine the total resistance of the conducting dia- 
grams. The remaining steps in calculating the diagrams 
are well known from the field theory of DP |2l], ^2| 23j . 

Renormalization group analysis provides us with the 
scaling behavior of the correlation function 



^ (O)^r(x) 



H 



il-d-r; 



■''II 



x jl + W A 2 M^/„ 
xfl + wX 2 x^f w ,, 



XL 



(6) 



where rj, v±, and z = vu/i/± are the critical exponents 
for DP known to second order in e-expansion [p2[ 2l| . fx, 
/2, fw,i, and f Wt 2 are scaling functions. In Eq. (|6|) we 
introduced the resistance exponent <j). Exploiting the 
fact that the correlation function in Eq. ([|) is a gen- 
erating function for Mr (cf. Eq. (||)) we deduce that 

Mr ~ |x|| l"^"" if measured parallel to the distinguished 
direction. For measurements in other directions it is ap- 
propriate to choose a length scale L and to express the 
longitudinal and the transverse coordinates in terms of 
L: |xjj ~ L and xn ~ L z . With this choice the scal- 
ing function f Wt x reduces to a constant and we obtain 
Mr ~ L^/ u± . For the resistance exponent we find in 
e-expansion 



= 1 



e 151 -314 In (4/3) 
24 + 6912 



+ O (e 3 ) 



(7) 



where e = 5 — d. Note that (f> is larger than the cor- 
responding resistance exponent for the random resistor 
network (RRN)|9| p5| . This is intuitively plausible since 
long tortuous paths that contribute to the macroscopic 
conductance in RRN are suppressed in DP. 

Now we compare Eq. (|7|) to the few numerical re- 
sults available in the literature. We are not aware of 
any numerical results for <j) itself. However, Redner 
and Mueller p6| determined the conductivity exponent 
t = <j) + (d — l)fj_ + Mi in two dimensions by Monte 
Carlo simulations: t(d = 2) = 0.6 ± 0.10. Arora et 
al. |27j did analogue and numerical simulations leading to 
t = 0.73 ± 0.10. Another value for comparison is t w 0.7 
obtained by Redner fl5]] from a real space renormalization 
group calculation. Crudely evaluating the e-expansion of 



3 



t for small spatial dimensions leads inevitably to poor 
quantitative predictions. Therefore it is appropriate to 
improve the expansion by incorporating rigorously known 
features. By carrying out a rational approximation which 
takes into account that t(d — 1) = we obtain the inter- 
polation formula 



t 



1 



(2 + 0.2083 e + 0.0604 e 2 



(8) 



which leads to t(d = 2) « 0.8. 

The second part of this letter is devoted to the back- 
bone dimension in DP. The backbone between two sites 
x and x' is defined, apart from unimportant Wheatstone- 
bridge-type configurations, as the union of all bonds car- 
rying current when / is inserted at x and withdrawn at 
x 1 . The average number of bonds belonging to the back- 
bone is referred to as its mass Mb- For the purpose of 
calculating the backbone dimension Db we assume that 
the bonds obey a generalized Ohm's law 



(Yj - Vi) [Vj - Vi\ [Vj - v t 



,1/r-l 



I, 



(9) 



The parameter r measures the nonlinearity of the bond. 
<7jj (Vj — Vi) takes on the same values as in the first 



part. The field theoretic Hamiltonian for the gener- 
alized RDN is given by Eq. (Eh with wX 2 replaced by 

In the following we are going to use that the average 
resistance of the generalized RDN, Mf> r , and the back- 
bone mass are related via Mb ~ lim r ^_ 1 + Affj r (see, e.g., 
Rcf. 11 ). A two-loop calculation analogous to that in 
part one reveals that the coupling proportional to w r does 
not require an individual renormalization in the limit 
r — > — 1 + . As a consequence we obtain 



lim (pr/vi 

r— > — 1 + 



V 



(10) 



at least to second order in e. From the generalized version 
of Eq. (||) we deduce that the backbone mass scales as 



M B = |xJ 



Ju 



(11) 



For self-affine objects the notion of fractal dimension 
is less straightforward than for self-similar objects. To 
determine the fractal dimension of the DP backbone one 
considers a d — 1-dimensional hyper-plane with an ori- 
entation perpendicular to xm. The cut is a self similar 
object with the fractal dimension 



dc.„t = D B - I 



(12) 



where Db is the local fractal dimension of the 
backbone [p8|. According to Eq. (^Tj) the mass of the cut 
scales like 



Me* = |x_ 



\4>-l/ v X _ — 1 



^|| fu 



(13) 



By choosing once more |xj_| ~ L and ~ L z we find 
that M cut ~ L<l>-x/vx-* n This leads via Eqs. @ and 
©to 



D B = l + 4>-i/vx-z = l-ri = d-2P/vx, (14) 

where j3 = vx{d — 1 + rj)/2 is the DP order parameter 
exponent known to second order in e|22], |23j| . 

Equation ( [l4] ) is in agreement with scaling 
arguments |2j| yielding that the fractal dimension 
of the transverse cut through a DP cluster with local 
dimension dj is df — 1 = d—l — /3/vx- The analogous cut 
through the backbone can be viewed as the intersection 
of the cut through the cluster and the clusters backward 
oriented pendant ]15|, 



27 



Hence, the codimension of 
the backbone cut is twice the codimension fijvx of the 
cluster cut, which leads again to Eq. (|l4|). 

We conclude with a few remarks. Our approach gives 
Eq. ( |l4| ) perturbatively to second order in e, while the 
scaling arguments leading to Eq. ( |l4] ) are exact. Hence, 
Eq. (nj) has a manifestation in the renormalization group 
framework in form of some Ward identity. The fact that 
w_i rcnormalizes trivially to two-loop order is reminis- 
cent of this Ward identity. It is an interesting issue for 
future work to identify the Ward identity and its under- 
lying symmetry. Our result for the resistance exponent <p 
is for dimensions close to five the most accurate analytic 
estimates that we know of. In two dimensions our results 
show reasonable agreement with the known numerical re- 
sults. It is certainly desirable to have more and firmer 
numerical data for comparison with our analytic results, 
in particular in three dimensions. We hope that this let- 
ter triggers further simulations of transport in DP. 
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